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' O ' Abstract. The long-time dynamics of the critical contact process which is brought 

^ I suddenly out of an uncorrclatcd initial state undergoes ageing in close analogy with 

(\ • quenched magnetic systems. In particular, we show through Monte Carlo simulations 

in one and two dimensions and through mean-field theory that time-translation 
*^ I invariance is broken and that dynamical scaling holds. We find that the autocorrelation 

^ ■ and autoresponse exponents Ar and Xr are equal but, in contrast to systems relaxing 

to equilibrium, the ageing exponents a and b are distinct. A recent proposal to define a 
non-equilibrium temperature through the short-time limit of the fluctuation-dissipation 
\^ . ratio is therefore not applicable. 

o 
o 
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H ' 1- Introduction 

>: 

^ . It is a well-established fact, theoretically known for many years and also experimentally 

j^ . observed, that the long-time properties of the ageing behaviour in glassy systems can be 

rationally interpreted in terms of a scaling picture which allows to make sense of many 
peculiarities of these systems which were once thought to be non-reproducible, see [T]. 
It has since been understood that similar effects also occur in non-glassy, e.g. simple 
ferromagnetic systems. From a microscopic point of view, the basic mechanism appears 
to be the formation of correlated domains of a time-dependent typical linear size i{t). 
For simple ferromagnets and other non-glassy systems, it is generally admitted that 
£{t) ~ t^/^, where z is the dynamical exponent. Furthermore, if 4>{t^r) stands for the 
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order parameter of such a simple ferromagnet, it turned out that these non-equihbrium 
relaxation phenomena are most conveniently studied through two-time quantities such 
as the two-time connected autocorrelator r(t, s) and the two-time linear autoresponse 
function R(t, s) defined by 



r(t,s) = (A0(t,r)A0(s,r)) , R{t,s) 



(1) 

h=0 



6h{s,r 

where A(f)(t,r) = (t>{t,r) — {4>(t)) are the time-dependent fluctuations of the order 
parameter and h is the magnetic field conjugate to 0. Causality implies that R{t, s) = 
for t < s. By definition, a system is said to undergo ageing, if T{t, s) or R{t, s) does not 
merely depend on the time difference t = t — s, but on both the observation time t and 
the waiting time s. For recent reviews, see El IH El IHl El • 

From experimental, numerical and analytical studies on ageing glasses and 
ferromagnets it has been learned that these systems may display dynamical scaling 
in the long-time hmit [HEl El IH El El El • Specifically, consider the two-time functions 
in the ageing regime t ^ t^icm, s ^ t^icro and t = t — s ^ ^micro, where tmicro is some 
microscopic time. Then one expects the scaling behaviour 

r(t, s) ~ s~'frit/s) , R{t, s) ~ s-'~^fn{t/s) (2) 

where the scaling functions fr,Riy) have the following asymptotic behaviour ioi y ^ oo 

My) ~ y-'^/^ , fniy) ~ y-'^/^ (3) 

Here Ar and Xr are called the autocorrelation jHl E] and autoresponse [12] exponents, 
respectively. 

Up to present, most systems studied were such that they relax towards an 
equilibrium steady-state. Here we wish to ask what aspects of the ageing phenomenology 
remain if time-dependent systems with a non-equilibrium steady-state are considered. 
Besides the obvious question concerning an eventual scaling behaviour of two-time 
observables, we shall be specifically interested in the following two points. 

(i) For systems with detailed balance, it is convenient to measure the distance of a 
system from equilibrium through the fluctuation- dissipation ratio ^^ 

X(t,.):=r/2(t,s)(^^) \ (4) 

At equilibrium, the fluctuation-dissipation theorem (FDT) states that X{t, s) = 1. 
Systems with an equilibrium steady-state may also be characterized through the 
limit fluctuation-dissipation ratio 

Xoo= lim (\imX{t,s)) . (5) 



-»oo \t~»oo 

If the ageing occurs at the critical temperature Tc, X^o should be an universal 
number fT^ IT^ E] and this has been confirmed in a large variety of systems in 
one and two space dimensions [iSl (111 ESI 1111 EI , whereas the value X^o = is 
expected for temperatures T < T^. The order of the limits in (j^l) is important, since 
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limt^oo {^i^^s->oo ^ {t , s)) = 1 always. On the other hand, Sastre et al. ^B] asserted 
recently that a genuinely non-equilibrium temperature might be defined through 



and confirmed this through explicit calculation in the 2D voter model as well as in 
other connected spin systems that do not satisfy the detailed balance condition. 
Does their definition (jU)) extend to more general non-equilibrium systems (in 
particular those exhibiting an absorbing phase)? 

(ii) For ageing simple magnets, it has been proposed that dynamical scaling might 
be generalized to the larger local scale-invariance ^H]- Such local scale- 
transformations, with dilatation factors depending locally on space and on time, 
may indeed be constructed for any given value of z. From the condition that 
i?(t, s) transforms covariantly under the action of local scale transformations it 

follows da cni 



l+a-Xu/z 



R{t,s)=rA-\ {t-s)-^-\ (7) 

where r^ is a normalization constant. This prediction has been confirmed in several 
models, notably the kinetic Ising model with Glauber dynamics, both in the bulk 
|19l 1201 ^^"i close to a free surface [25 , the kinetic XY model with a non-conserved 
order parameter |^ 1^ . for the Hilhorst-van Leeuwen model J2I] and several 
variants of the exactly solvable spherical model jTHl CSl 1211 122] •§ Very recently, 
the autocorrelation r(t, s) could be predicted [221 from local scale-invariance for 
phase-ordering, where z = 2. This prediction has been confirmed in the Glauber- 
Ising, kinetic spherical and critical voter models and also for the free random walk 
123 122] • For recent reviews, see pHjIH]. 

Is there further evidence in favour of local scale-invariance in statistical systems 
without detailed balance ? 

Probably the simplest kinetic system far from equilibrium and without detailed balance 
is the celebrated contact process which has a steady-state transition in the directed 
percolation universality class. We shall therefore use this model in order to gain insight 
into the role of conditions such as detailed balance into the phenomenology of ageing 
behaviour. A complementary paper studies the same model through the density-matrix 
renormalization group [23 • In section 2, we recall the definition of the model. In section 
3, we discuss the computation of correlators either for the non-critical system or at the 
critical point. In section 4, we define the response function and estimate it at criticality. 
After that, in section 5, we discuss the possibility of a generalization of the FDT for this 
model and in particular whether a non-equilibrium temperature may be defined through 
eq. ©. Finally, we conclude in section 6. In the appendix, the mean-field theory of 
ageing in the contact process is discussed. 

§ A technical complication arises for the \D Glaubcr-Ising model where the construction of the Lie- 
algebra generators of local scale-invariance must be generalized, see appendix C in |23j . 
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2. The model 

The contact process (CP) was originaUy conceived as a simple model to describe 
epidemic disease propagation j2Hl 121] • Since then, CP has become a paradigm 
of Directed Percolation (DP), arguably the most common universality class of 
nonequilibrium phase transitions with an absorbing state. It is precisely the existence 
of an absorbing state that ensures that this model does not satisfy the detailed balance 
condition. This property stimulates our interest to try and characterize the model's 
dynamics in a new way. 

We begin by recalling the definition of the model. In the contact process, the states 
of the system are described by a discrete variable, rij(t), defined on the sites 2 of a 
hypercubic lattice. The possible values of rzj are 1 or depending on whether the site i 
is occupied or empty. The dynamics is defined as follows: for each time step, a site i of 
the lattice is randomly selected. If i is occupied, that particle vanishes with probability 
p. Otherwise, with probability 1 — p, a. new particle is created on one of the nearest 
neighbours of i chosen at random (if that new site was still empty) . When the control 
parameter p is varied, the model exhibits a continuous phase transition from an active 
phase, where the mean density {n{t)) tends to a constant value n in the stationary state, 
to an absorbing phase, with zero final density. Separating these two phases, there is 
a critical point located at pc = 0.2326746(5) in ID jHOl IHIj (the number in brackets 
gives the uncertainty in the last digit), and at Pc = 0.37753(1) in 2D jSI]. We have 
arbitrarily fixed the initial condition of our system such that the particles are randomly 
distributed throughout the lattice with a mean initial density of Uq = {n{0)) = 0.8. 
The linear system sizes used are L = 10^ in ID and L = 300 in 2D, respectively; the 
typical number of disorder realizations considered in simulations goes up to 2 x 10^ in 
one dimension and 7 x 10^ in 2D. 

The introduction of an external field in the model will be necessary to estimate 
a response function. For this model, the external field h corresponds to an additional 
probability of creation or destruction of particles, depending on the sign of h. 

3. Two-time autocorrelation functions 

A crucial difference between the CP and usual spin models is that starting from a zero- 
density initial state (corresponding to zero magnetization in magnets) is not possible. 
This configuration corresponds to the absorbing state of the CP and once the system 
falls there, no further evolution occurs. Therefore, the spatially averaged mean density 
{n(t)) is positive during the time the systems evolves. For example, it decays as a power 
law, {n(t)) ~ t~^, at the critical point p = Pc- This fact implies that the two-time 
autocorrelation function defined as 

C{t,s) = {n{t,r)n{s,r)), (8) 
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Figure 1. Connected autocorrelation function r(t, s) for the CP for several values of 
the waiting time s. In a) the steady-state is in the absorbing phase. The main plot 
shows the ID system with p ~ 0.286 and the inset shows the 2D system with p = 0.56. 
In b) the steady-state is in the active phase with p = 0.2 in the main plot for ID and 
p = 0.286 in the inset for 2D. The system sizes are L = 10^ in ID and L = 300 in 
two dimensions. The waiting times in a) are from bottom to top s ~ [30, 100, 300] and 
s ~ [3, 10, 30] in the inset. In b), they are, from bottom to top, s ~ [1, 10, 30, 100, 1000] 
and s = [0.1, 1, 3, 10, 100] in the inset. 



which is widely used in spin systems or glassy models, is no longer equal to the connected 
autocorrelation function or covariance that reads 

T{t,s) = {An{t,r)Anis,r)), (9) 

where An{t, r) = n{t, r) — {n{t)). Both functions are related by means of the expression 

T{t,s) = C{t,s)-{n{t)) {n{s)). (10) 

In eq. fllOj) . the second term on the right behaves asymptotically as 



{n{t)) {n{s)) 



(ts) 



iip = Pc 
Up < Pc 



:iii 



Therefore, for p < Pcit follows that C{t, s) and T{t, s) do not share the same asymptotic 
behaviour. 

Once the correlations have been defined, we are in position to deal with a first 
physical aspect of the temporal evolution of the contact process. It is well-known that 
away from the critical point, the model has a finite relaxation time towards its steady- 
state. That is, the mean density decays exponentially fast to its stationary value n: zero 
in the absorbing phase or a positive quantity in the active one, see |2H]- In Figure^, the 
transient behaviour of r(t, s) for t — s relatively small is shown for a system evolving to 
the steady-state in the absorbing phase. We see that, not only the density, but also its 
fluctuations with respect to (n) are vanishing, as described by the drop of the F-curves. 
The signature of a finite relaxation time may be noted in the persistence of the curve 
shape throughout the system evolution. || Next, we consider whether a similar behaviour 

II Far above Pc, one has for large s C{t, s) ^ r(i, s) ^ {n{t)), see [2] for details. 
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Figure 2. Microscopic evolution of clusters in the critical 2D contact process. The 
lattice size is, in both series, 1000 x 1000. The initial condition of the upper plots is a full 
circle with radius 100 placed in the center of the lattice, while below it is a full lattice. 
The times are t = [2, 20, 200, 2000] for the upper figures and t = [20, 200, 2000, 20000] 
for the bottom snapshots. 



can be seen in the active phase as well. In Figure^, the connected correlation function 
is displayed for a system in the active phase in both ID and 2D. The curves move 
upward during a first period of waiting times before they collapse, which implies that 
the TTI condition is fulfilled in this phase. 

We point out that the observation of TTI at large enough waiting times in the 
active phase of the contact process is qualitatively different from what is found in 
ferromagnets quenched into their ordered phase. This may be understood qualitatively 
from the Ginzburg-Landau functional by observing that for ferromagnets there are two 
stable competing steady-states while for the contact process there is only a single one, 
to which the system relaxes rapidly. 

From now on, we concentrate on the behaviour at criticality. In spin systems with a 
global symmetry (e.g. an up-down symmetry in the Ising model), the critical dynamics 
is based on the growth of spin domains whose average size, associated to the correlation 
length, increases as a power law in time, i{t) ~ t^/^ where z is the dynamical exponent. 
The time-evolution of typical configurations of the critical CP in two dimensions is 
shown in Figure |21 Note the substantial difference with respect to bulk ferromagnetic 
spin systems: no domain walls are formed and, at the same time, the number of occupied 
sites decreases in time. This kind of ageing behaviour, not driven by a surface tension, 
has already been described in the 2D voter model universality class ^7\. Recently, it 
was pointed out that cluster dilution may also occur in the early stages of surface ageing 
of magnetic systems |2H]- As we shall see, the asymptotic dilution of clusters leads to 
new features in the scaling description of the ageing process of the CP. 

In analogy with the behaviour of magnetic systems at the critical point, a scaling 
ansatz of the type 
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T{t,s) 
Cit,s) 



'fcit/s) , fc{y)-^y~^^'' ify^oo 



(12a) 
(12b) 



may be expected for the correlation functions, where the exponents 6, \c/ z and \y/ z 
are to be determined. In Figure Et^, the critical evolution of r(t, s) is plotted for several 
values of the waiting time s, over against t — s. Time-translation invariance is now 
clearly broken, and we find ageing effects for both one- and two-dimensional systems. 
The data can be collapsed when plotted over against t/s, as shown in Figure Eb which 
provides evidence for the validity of the scaling relation eq. p2h .). The values of the 
exponent h used for the collapse were obtained from the scaling relation 

(rf - 2 + r/) 






2 5, 



(13) 



which remains valid for the DP universality class. Here, j3 is the static exponent which 
describes the variation of the order parameter close to criticality n ~ (pc — p)^, ^|| is the 
temporal correlation-length exponent, r] is the static exponent of the order parameter 
spatial correlations, 6 is the exponent controlling the critical decrease of the mean density 
and z is the dynamic exponent (their values are summarized in table 1). 

The exponent Xr/z may be estimated from the collapse of Figure Eb? with 
reasonable numerical accuracy. We find Xr/z = 1.9±0.1 in ID and Xr/z = 2.8±0.3 
in 2D. If the critical CP were a Markov process, these exponents might be calculated 
from the global persistence exponent 6g by means of the scaling relation jSHl EDI 1^ 



Ar „ 2{l-d)~v 

z 2z 



(14) 



This expression predicts Xr/z = 1.98(2) in ID and Xr/z = 3.5(5) in 2D, not too far 
from our numerical measurements. Still, the deviations from the directly measured 
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Figure 3. a) Evolution of the connected autocorrelation function at criticality. The 
curves in the inset are for two dimensions and the main plots for ID. The waiting 
times are in both cases s ^ [10, 20, 50, 100, 200, 500, 1000, 2000] from left to right. In 
b) the collapse of the previous curves using Eqs. H12h , ll3() is shown. The straight lines 
yield exponent values of \r/z ~ 1-9 in ID and of Xr/z = 2.8 in 2D. 
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Figure 4. a) Two-time autocorrelation function C{t, s) for the contact process in 
ID (main plots) and in 2D (insets). The waiting times arc from top to bottom 
s = [100, 200, 500, 1000] in both cases; for one dimension the value s ~ 2000 is also 
included. These curves are collapsed in b) employing the scaling ansatz of Eq. H12bll3|l . 
The straight lines have a slope —6. 



values appear to be significant, pointing towards the possible existence of temporal 
long-range correlations and, hence, of an effective non-markovian dynamics. 

The evolution of C{t, s) is shown in Figure ID for a range of waiting times. The 
collapse of these curves with the scaling ansatz of Eq. ()12p is also shown. From this 
collapse, we make the identification Xc/z = 6 and also confirm the relation b = 26. As 
expected, we have found a different asymptotic behaviour for T(t, s) and C{t, s). These 
results in ID and 2D can be compared with the mean-field prediction in the scaling 
limit s ^ oo with y = t/s kept fixed (see the appendix) 



Cit,s)^- 



(15) 



This is in perfect agreement with the scaling ansatz eq. (|12b) since the mean-field value 
for the exponent S is 6mf = 1 (see table 1). 



4. The two-time autoresponse function 

The linear autoresponse function at time t to an external perturbation by the field h at 
time s is defined as 

6{nit)) 



R{t,s) 



Sh{s) 



(16) 



h=0 



Unfortunately, this function is very difficult to estimate directly. Hence, as a general 
rule, the susceptibility is measured instead. It is common to work with the zero-field- 
cooled (ZFC) susceptibility (we keep the terminology of magnetic spin systems by abuse 
of language) 

XzFc{t,s) =x{t,s) = / duR{t,u). (17) 
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Table 1. Exponents for the directed percolation universality class in one and two 
dimensions and in mean-field theory. The exponent b = 26 and furthermore Xc /z = S. 
The numbers in brackets give the estimated uncertainty in the last digit (s). 



To measure tliis, one may consider tlie following dynamical rules. First, the system 
evolves according to the original CP rules, until the waiting time s has elapsed. After 
that, two copies of the systems, A and B, are kept evolving in parallel. The copy A 
continues its evolution without perturbation, while copy B is subjected to an external 
field h. This field, as explained in section 2, represents a certain probability of creation 
{h > 0) or destruction {h < 0) of a particle whenever the selected site is empty or 
full, respectively. The susceptibility is then calculated as the limit h —>■ oi the ratio 
between the difference of the densities of the two copies and h: 

where h is turned on at time s and kept until time t. However, the application of this 
method to the contact process is not straightforward, since a field h ^ suppresses the 
phase-transition and hence the ageing behaviour. For spin systems, this difficulty may 
be side-stepped by using a spatially random magnetic field hi = ±h [S21- Although 
this method may be applied to CP, the random field still brings the system out of its 
critical point, since there is no symmetry between the two states of the variable Ui. 
Consequently, to measure the susceptibility, there is no option but to force the system 
out of criticality. We have checked that applying an uniform field to the whole system 
is numerically more efficient than utilizing a random one. Therefore we shall employ 
the former method in our calculations. We stress that, in both procedures, a small 
enough value of h must be taken to prevent saturation for the times considered in the 
simulations. From studies of the magnetization-reversal transition in spin systems it 
is known that ageing behaviour can be found if h is small enough |221- In Figure Et, 
the susceptibility is depicted for a fixed waiting time and several values of the external 
field. The asymptotic behaviour {h — > 0), which corresponds to the early r = t — s 
regime, may be there clearly observed. In addition, in Fig. Eb, the same magnitude 
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Figure 5. a) The ZFC susceptibility x{t, s) as a function of t — s for the one- 
dimensional contact process ai p ~ pc- The external field is changed for the different 
curves, from bottom to top h = [10^^, lO^'', 10^"*, 10"^], while the waiting time is 
maintained at s = 1000. In b), the external field is fixed at h = 10""^ and x(^7 5) is 
represented for several waiting times, from bottom to top, s = [250, 500, 1000, 2000]. 



is represented for a range of waiting times. This shows that there is a time window 
available for studies of the ageing behaviour of the linear response. 

The interpretation of the zero-field-cooled susceptibility in terms of the scaling 
behaviour of the linear autoresponse function R{t, s) meets with a further difficulty. 
Naively, one would simply insert the assumed scaling form eq. into (|T7jl . This 
simplistic procedure would lead to 

X{t,s)= [ duu-'-yR{t/u) = t-^ [ \ww''-^fR{w) = s-''fM{t/s) (19) 



However, as pointed out in [Mj, in doing so one neglects the important condition 
t — s ^ ^micro necessary for the validity of eq. ^. Indeed, it can be shown that taking 
this condition into account rather leads to j31] 



X(t, s) = xo + s~^gM{t/s) + O (s"*^) 



(20) 



such that the exponent A > is in general unrelated to a and where Xo is some constant. 
For example, for phase-ordering spin systems which are at a temperature above their 
roughening temperature Tr, one always has A < a j^l]. Hence the term of order 0(s~°) 
in pn|l . coming from the integral (|TTI|l over the autoresponse function, merely furnishes 
a finite-time correction (see the appendix for a related difficulty of using Xzfc which 
arises in mean- field theory). 

These difficulties may be circumvented by using the integrated response in a 
setup similar to a Thermoremanent Magnetization (TRM) experiment, with a small 
modification, close in spirit to the 'intermediate' protocol proposed in [SI]. To do so, 
we start with a system evolving with the usual contact process rules. After a time 
s — Tc, the system is split in two copies, A and B. A continues its evolution without 
perturbation, whereas B is subjected to an external field h. The field is switched off at 
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time s, and from then on we track the difference between the two copies, which is called 
the thermoremanent density: 

Pit, s) = r du R{t, u) = hm Mt)-nAit)) ^ ^21) 

/„ ^ h^O \h\ 

'J S — Tc II 

where h is now turned on at time s — r^ and switched off at time s. In addition, t^ must 
satisfy the relation Tc< s - actually, this scheme would correspond to a standard TRM 
experiment if Tc = s. In addition, Tc must also be much smaller than the relaxation time 
of copy B to avoid that p(t, s) becomes a function of the difference of times t — s alone 
j32] • The scaling of p may be obtained from inserting eq. ((Tj) of local scale- invariance into 
fl21|) . It can be shown [S^ that the leading corrections to scaling are of order 0(5"^^/^) 
which is negligible at criticality. 

A last point has to be discussed before the scaling of the integrated response p as 
calculated here can be analysed quantitatively. Consider the linear space-time response 
function 

5{n{t,r)) 



R = R{t, s;r- r') 



(22) 

h=0 



5h{s,r') 

In contrast to the situation usually considered in magnetism, where the field h is 
modelled by spatially uncorrelated random variables, we use here an uniform magnetic 
field with spatial correlation, as explained before. This implies that the measured TRM 
becomes for s ^ Tc 

P(t^ ^) = ^J2Y1 ^(^' "' ^ - ^') = ^o(^' ")' (23) 

r r' 

where the sums runs over the entire lattice of linear size L and Rk{t, s) is the Fourier 
transform of R{t, s; r) at momentum k. In order to understand the scaling behaviour 
of Ro(t,s), we recall that local scale-invariance predicts for any given value of z the 
following form of the space-time response ^H] 

R{t,s-r) = R{t,s)'^{\r\/0, (24) 

where R{t,s) is the autoresponse function, ^ ~ (t — sY^^ and $ is a scaling function 
which can be obtained from a linear differential equation of fractional order. ^ Therefore 

Roit, s)= f dr R{t, s; r) ~ R{t, s)^'^ ~ R{t, s){t - s)'^/^ (25) 

where we assumed that the scaling function $(m) falls off sufficiently rapidly for u 
large such that the Fourier transform exists. If the autoresponse function scales 
asymptotically (i.e. for t/s ^ 1) as 

i?(t,.)~s-i-M- (26) 



^s ^ 
we expect for the TRM the scaling behaviour 

p(t,s)~^o(t,s)~s-i"'^+'^/M- (27) 
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Figure 6. a) Temporal evolution of the 'thermoremanent' density p{t, s) (see eq. (I21() ') 
for several waiting times, from bottom to top s = [250, 500, 1000, 2000]. The main plots 
are for one dimensional systems and the insets for 2D. In b), we collapse the curves of 
a) according to dynamical scaling cq. (|27|l . The slopes of the straight lines are —1.27 
in ID and -1.62 in 2D. 



In Figure ^l, the evolution of p is plotted for several values of the waiting time. 
The strength of the field used in the simulation is h = 10~^ and Tc = 25. The ageing 
behaviour is clearly visible. Replotting the data according to the dynamical scaling 
behaviour expected from eq. (jTfjl (figure IHb), we find a nice data collapse for the 
exponent values 



1 + a 



d _ I -0.20(10) ; in ID 
i ~ j 0.17(10) ; in 2D 

Taking the slopes of the curves in figure IHb we get 

XR-d I 1.27(10) ; in ID 
1.62(10) ; in 2D 



(28) 



(29) 



from which the exponent values of a and Xr given in table 1 are obtained. 

It is satisfying that for d = 1, the results of a complementary DMRG calculation 
fI7\ are consistent with ours. Our results may also be compared with the prediction of 
mean-field theory, as derived in the appendix 

-2 



R{t,s) 



it 



~d/2 



(30) 



The mean-field values for the exponents are then a = | — 1 and Xr/z = | + 2. They 
should become exact ior d > 4 and are not too far from the simulational results in lower 
dimensions. 

Comparing with the exponents obtained from the autocorrelation function in 
section 3, we observe that, within numerical errors, the autocorrelation and autoresponse 
exponents agree 

Ar = Xr, (31) 

% For z = 2, it can be shown that $(u) = exp(— A^it^). where A1 is a dimensionful constant |18| . 
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but that the exponents a and b are different. In fact, within error bars they obey the 
following relation 



26. 



(32) 



This last result leads to important consequences as we shall see in the next section. 
Finally, the uncertainty in the ageing exponents derived from our numerical data is 
too large for a direct check of the validity of local scale-invariance to be carried out. 
Nevertheless, our results can be related to the ones using a direct estimation of R{t, s) via 
the DMRG J2ZI through the ansatz ()24j) . which is based on local scale-invariance. The 
agreement between the two measurements provides indirect evidence for the existence of 
this symmetry in the CP. On the other hand, the prediction ((Zj) of local scale-invariance 
is satisfied by our mean-field result (fHOjl and in ID, the vahdity of eq. ((Zj) was directly 
confirmed from a DMRG calculation fI7\ . 

5. Fluctuation-dissipation theorem 

Our observation that a 7^ fe at criticality has important consequences. In particular, it 
follows that no analogue of the fluctuation-dissipation ratio can be defined in the critical 
contact process in a meaningful way. To see this, we use the scaling ansatz (fT^ ) for 
the connected correlation function and find 
dT{t,s) 



ds 



-'[bMt/s) + {t/s)f^{t/s)]. 



(33) 



where the prime denotes the derivative. This implies for a formally defined fluctuation- 
dissipation ratio R(t,s){dr(t,s)/ds)~^: 

(i) a limit fluctuation-dissipation ratio, in analogy with (jSJ, can only have a non-trivial 
value if simultaneously a = b and Xr = Ap. From table 1, we see that this is not 
the case for the critical contact process. 



200 



CO 




0.1 0.2 

r(t,s) 



Figure 7. The ZFC susceptibility against the connected correlation function for 
a one-dimensional system and with different waiting times, from bottom to top 
s = [250, 500, 1000, 2000, 5000]. The external field was fixed at /i = 7 x 10"^. 
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(ii) in contrast to the finding of Sastre et al. jTHI for the critical voter model, there is 
no non-trivial limit t — s^O inside the scaling regime of the fluctuation-dissipation 
ratio. Existence of such a limit requires a = b a.s a. necessary condition. From table 
1, that condition is not satisfied. Consequently, the supposedly generic definition of 
a non-equihbrium temperature -eq. ^- as proposed in (TH] does not extend to the 
critical contact process and hence should reflect a peculiarity of the models with 
symmetric states."*" 

A more graphical way to arrive at the same conclusion is to plot the ZFC 
susceptibility xit^^) versus the connected correlation function r(t, s): if there were 
a generalized FDT, we should find a curve with a constant slope that corresponds to 
the inverse of the effective temperature. We do so in figure [7| for the ID case, which 
makes it clear that no such extension of the FDT exists, in contrast with results fH] of 
the 2D critical voter model. 

6. Conclusions 

We have explored the role of the detailed balance condition in the slow non-equilibrium 
relaxation of a statistical system towards a steady-state. In order to do so, we have 
replaced a key property of physical ageing, namely that the basic microscopic processes 
are reversible, by a property encountered in chemical or biological ageing with their 
underlying irreversible microscopic processes. Consequently, the steady-state of the 
system can no longer be at equilibrium. One of our main questions was how the 
accepted scaling behaviour found in physical ageing might be modified in the absence of 
detailed balance and we tried to get some insight into this through a case study of the 
contact process, which is a paradigmatic example of systems with a steady-state phase- 
transition in the directed percolation universality class. We have studied the ageing of 
this model by analyzing two-time autocorrelation and response functions, obtained from 
Monte Carlo simulations and also analytically from the mean-field approximation. Our 
numerical results are in full agreement with those coming from an application of the 
density-matrix renormalization group to the same model |27] . 

Our results reemphasize the importance of having either at least two physically 
distinct and stable steady-states or else being at a continuous phase-transition in order 
to be able to observe the slow non-exponential relaxation towards a global steady-state 
which is associated with the phenomenological aspects of ageing behaviour, notably the 
breaking of time-translation invariance and of dynamical scaling. In this sense, it is 
natural that in the contact process ageing is only found if the steady-state is critical. 
We have shown that the phenomenological scaling description, previously developed 
for ageing magnetic systems, generalizes to the case at hand. In particular, we have 

+ The basic of assumption of ^HIj namely that X{t, s) ^ 1 in the short-time hmit has been criticaUy 
reexamined by Mayer and Solhch 0SI, who in the ID Glauber-Ising model at T = construct a 
defect-pair observable such that lim^^^oo X{s, s) = 3/4. 
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presented evidence for the exponent relations 

Ar = Aij , l + a = b = 26 (34) 

in one and two dimensions and in mean-field theory. While the first of those is generic 
in ageing magnets with an uncorrelated initial state, the second one is not and suggests 
that the non- equilibrium exponents a andb are in general distinct* It will be interesting 
to see whether or not the relations ()34|1 hold at the critical points of other nonequilibrium 
universality classes. Eq. (pij) implies that a recent attempt ^H] fo define a non- 
equilibrium temperature through an extension of the fluctuation-dissipation theorem 
cannot be generalized beyond the rather restricted context where this definition had 
been proposed. 
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Appendix. Mean-field theory 

The contact process is described at the continuum level by the Reggeon Field Theory 
(RFT) j33]. Let (f){t, r) represent the continuum order-parameter at the spatial position 
r and at time t. Taking into account the initial conditions, which is fundamental to 
analyze this problem, the RFT-action reads 

S = fdtdr {dt(p -met)- DV^cf)) + gcj)"^^ - gcl)^)'^ - noS{t)4> , (Al) 

where is an auxiliary field, m the mass term, D the diffusion constant, g a coupling 
constant and uq the initial density. The mean-field approximation amounts to treating 
2]) as a classical action. This yields the following equations of motion 
6S 



S<Pit) 

and 

5S 



{dt-m- DV^) + ^0^ - no5{t) = (A2) 



{-dt -m- DV^) + 2^00 = 0. (A3) 



We can now evaluate the classical density and the classical response function. The term 
classical means here averaged with respect to action (jAljl but only with the 0^0 vertex 
(diagrammatically, the term classical refers to the absence of loops in the diagrams). 
The classical density is given by the sum of all tree-diagrams which terminate with a 

* For systems with detailed balance, the steady-state is an equilibrium state where the fluctuation- 
dissipation theorem holds. This in turn implies a = b |27| . 
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single propagator. If one evaluates these diagrams in momentum space, it follows that 
the 6{t)(j){t) only make a contribution with k = and therefore all diagrams at tree level 
have A; = 0, see jlH]. The classical density is found by Fourier-transforming eq. ()A2|1 . 
We find 

k^) = — ^ + ^- / ^k^ - ^')k^% (A4) 

m + luj m + 1UJ J Zn 

where 4>{oo) is the Fourier transform of (f){t). Note that the diffusion term has been 
disregarded since it gives no contribution to the final result, as explained above. Eq. 
()A4|) can be iterated to obtain a perturbative series. For m = 0, that is, at the critical 
point, the series can be summed yielding the exact solution 

m = -r-^ (A5) 

1 + gnot 

Therefore, the classical two-time correlation function is given by 

C{t, s) = 0(t)0(s) = r^- -. (A6) 

We can now turn to the computation of the response function. The Fourier 
transform of equation ()A3|1 with respect to r reads 

- (9t0(fc, t) = m4>{k, t) - Dk^4>{k, t) - 2^0(t)0(fc, t), (A7) 

where 4'{k,t) is the Fourier transform of (j){t,r). The above equation has the formal 
solution 

0(fc,t) = 0(fc,s)e(-"^+^'=')(*-^)exp 2g f dt' (f){t') . (A8) 

•J s 

At criticality, m = and using equation ()A5j) the above equation yields (setting D = 1) 

0(fc,t) = 0(fc,s)e'»^(-)fl±^y. (A9) 

\l + gnosJ 

In order to calculate the response function is necessary to take t > s and to multiply 
both sides of equation ()A9|) by (f){—k,t + e), where e > 0. The fact that there is no 
inertial term in the action (jAlj) (that is, there are no second-order time derivatives) 
leads to the following boundary condition 

lim /0(-fc, t + e)0(fc, t)\ = \imRk{t + e, t) = 1 (AlO) 

Consequently, the classical response function reads 

(-fc, t)0(fc, s)) = R,{t, s) = e-^^(*-) (^^^) ' (All) 

In real space, the space-time-dependent response function becomes 

R{t, s; r) = (0(t, r)0(s, 0)) = i?(t, s) exp ('-—!—-') (A12) 

where the autoresponse function is 
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We are interested in these expressions in the ageing regime s^l,t — s^l. The 
the autocorrelation function reads 

C{t,s)r^s~H-' = s~^(-\ (A14) 

and we read off the mean-field values 6mf = 2 and {Xc/^)mf = 2 of the exponents 
(mean- field theory being gaussian, zuf = 2 is evident). The autoresponse function 
becomes 

i?(t,s) = 7^'^/2(t-s)"'^/2|'i^ ^A^5) 

in agreement with the form predicted by local scale-invariance and we have the exponent 
omf = f ~ 1 and (Aj^;/2;)mf = f + 2. While simple mean-field does not yield a prediction 
for the connected correlation function r(t, s), at least it furnishes expressions for the 
correlator and the linear response and their exponents which should become exact for 
c/>4. 

We finally illustrate the importance of taking both conditions s ^ 1 and t — s ^ 1, 
or equivalently s ^ 1 and t/s ^ 1, for the validity of the scaling form ((2)) into account. 
First, for the thermoremanent protocol (we take Tc = s) we find, using ()A15j) 

p^^^{t,s)= duR{t,u)^t^-'^/^ Aww^l-wy^'^ (A16) 

which produces a well-defined scaling function for all values of d. On the other hand, 
for the zero-field-cooled protocol we have 

XzFc(t, s)= f du R{t, u) ~ t^-'^/' / dw w\l - wy/^ (A17) 

Js Js/t 

which diverges for all rf > 2. This means that the scaling form ()A15|1 cannot be used near 
the upper integration limit in ()A17j) and the contribution of the region around w = 1 
must be evaluated separately, without appealing to dynamical scaling (we remind the 
reader that the condition t — s ^ 1 is implicit in the continuum limit used to derive 

dnu). 

After submission of this work, we became aware of a paper by Oerding and van 
Wijland [IHI where the two-time connected autocorrelation function ro(t, s) of reggeon 
field-theory is worked out in momentum space to one-loop order. They find a scaling 
behaviour quite analogous to the one observed here. 
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